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I report some progress that occurred since NO-VE 08 in the field of non-standard 
d \ interactions (NSI) of neutrinos. After briefly reviewing theoretical developments, 

I give a summary of the two works in which I was involved. Firstly, we have 
formulated a perturbative framework to illuminate the global features of neu- 
trino oscillations with NSI, aiming at exploring method for determination of the 
standard mixing and the NSI parameters. We have recognized that the parame- 
f~i ■ ter degeneracy prevails with an extended form which involves the NSI elements. 

Oh| Furthermore, a completely new type of degeneracy is shown to exist. The nature 

^ • of the former degeneracy is analyzed in detail in the second work. The work is 

O \ primarily devoted to analyze the problem of discriminating the two CP violation, 

i-C i one due to the lepton Kobayashi-Maskawa phase and the other by phase of 

the NSI elements. We have shown that the near (3000 km)— far (7000 km) two 
detector setting in neutrino factory does have the discrimination capability and 
^ ' is sensitivities to CP violation due to NSI to \£efi\ to 
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of the region of 5 and (j)^^ ■ 
1. Introduction 



in 

^ . The question I would like to address in my talk is: Are there something terribly 

new in neutrino properties after the discovery of neutrino masses and lepton mixing? 
The Japanese saying in the subtitle is meant to be thatQ Clearly this is an extremely 
interesting question. But, I must start with a cautionary remark. 



I What is the natural time scale for discovery of something extremely new in neu- 

trino properties? Let us look back the history to obtain a hint for answering the 
question. It took more than 60 years from the Meitner-Hahn measurement of elec- 
tron energy spectrum in nuclear beta decay in 1911 to the discovery of NC reaction 
in 1973 From the discovery of neutrino itself in 1953 by Reines and Cowan ^ 
to the discovery of neutrino mass and lepton flavor mixing by Super-Kamiokande in 
1998® needed 45 years. Thus, the right time scale, as history tells us, is ~ 50 years. 
It implies a warning; What people think about the possible candidates for "terribly 



^ Written version of a talk presented at XIII International Workshop on Neutrino Telescopes, 
Istituto Veneto di Scienze, Lettere ed Arti, Venice, Italy, March 10-13, 2009. 

^ Though the original saying in Japan is "Another loach under a willow?" I decided to keep the 
modified version because eel is much more familiar to us, and in particular thanks to a great support 
to the eel version by the spokesman of IceCube collaboration. 



new in the neutrino properties" at the right time can be very different from those we 
consider today. Though I have to speak within the scope that I can think of today 
but this point has to be kept in mind as a word of caution. 

It also has to be remarked that the scope of my presentation is very hmited; 
Thou gh I res trict myself into the so called non-standard interactions (NSI) of neu- 
trinos l ^ ^ ^ l ^ l ^ l' in my talk, "new in the neutrino properties" may include various often 
more radical possibilities such as: departure from three flavor mixing, sterile neutri- 
nos, violation of fundamental symmetries like CPT. See e.g.,^ for a status summary 
for these more exotic possibilities. 



2. Non-Standard Interactions of Neutrinos 



It has b een pro posed that neutrinos might possess yet unknown new neutrino 
interactions I ^ I S I T lHI Today we have Standard Model of particle physics, one of the 
most successful theories in physics. Therefore, when we discuss NSI it is natural (and 
to a large extent mandatory) to talk about it in a language of higher dimensional 
operators 1^^. 

Suppose that there exist a new physics at energy scale M^p, which I assume to be 
greater than ~ 1 TeV, but not too much larger than this value. I assume the ty pe of 
higher-dimentional operators for effective interactions of neutrinos with matter EUD 

£NgSI = _2v^egG^(i7„7^P^z.^) (fYPf), (1) 

where Gp is the Fermi constant, and / stands for the index running over fermion 
species in the earth, / = e,u, d, in which we follow the conventional notation. P 
stands for a projection operator and is either Pp = 1(1 — 75) or Pp = ^(1-1-75). Given 
the dimension six operator in ([T]) and because we normalize the operator with Fermi 
constant Gp, 8^/3 must be of the order of (Mw/Mnp)^ ~ 0.01 (10"^) if M^p = 1(10) 
TeV. If we have to go to dimension eight operators their effective strength would be 
at most (Mw/Mnp)^ ~ 10"^ even for Mnp = 1 TeV. The off-diagonal elements may 
have further suppression. 

Since I gave a talk on NSI last year in Venice I will restrict myself into 
developments that occurred after NO-VE 08 to show that the field is moving. The 
rest of my report has three parts: In section 3 I review the recent development in the 
theory of NSI. From section 4 I change gear to NSI effect in propagation in matter. In 
section 5 I discuss perturbative treatment of the system with NSI. Sections 6 and 7 are 
devoted to further clarifying the properties of the system and to discuss the question 
of discriminating two kind of CP violation, one from the lepton Kobayashi-Maskawa 
phase and the other from phases of the NSI elements, the problem discussed in 
^1^. My presentation in the last three sections will be based on the two recent papers 



I ^^ l ^^ ^ The latter work is a natural continuation of our previous work^^^. 

3. Recent Development in the Theory of NSI 

Let me start by reviewing the development in the theory of NSI that occurred 
very recently. Since long time ago, it has been noticed that phenomenological 
study with NSI of the type ([T]) has a potential caveat. To get to the point, let us 
agree on the following understanding: At a high-energy scale where NSI originates 
the SU{2) X f/(l) gauge invariance holds. Then, the left-handed neutrino field in the 
operator ([1]) must be elevated into the lepton doublet of SU{2). When we require 
this an obvious problem occurs; The resultant four charged lepton operators have to 
obtain severe constraints from experiments. The most stringent is the one imposed 
by the branching ratio of /i — eee, BR{fi eee) < 10^^^ , which would yield 
the constraint \Ben\ < 10^^. Of course, nothing is wrong with it. But, we would 
like to avoid this because we are interested in observable effects in near (or even 
remote) future neutrino experiments. Therefore, people looked for the possible higher 
dimensional operators which are free from the charged lepton constrai nts. 

Some candidates which were discussed by people (see, for example, I 19 | 20 | 21 t^ g^^^. 

Ol = {L.iT^LDmiT^Ls) (2) 

for dimension six operator where L'^ = CL^ and C is the charge conjugation operator. 
For dimension eight operators they are of the type 

O^s = {Lp-ipL^){LsH)Y{H^L,). (3) 

See, for example, | 11 | 22 | 21|23|24 | relevant references. Intuitive understanding of ([3]) 
is that the Higgs field v.e.v. projects out only the neutrino component of left-handed 
doublet. The meaning of ([2]) becomes clear by writing it in a form with obvious 
antisymmetry in flavor space , 

201 = (47%)('^77M^5) + (^V^5)('^a7M^,3)-(47%)('>77M^/3)-(^"77'^^,3)(^>a7M^5) (4) 

which implies 

a/3 _ 7/3 _ a5 _ 7'5 /cN 

The antisymmetric nature prohibits, for example, e^^^ which would produce NSI effects 
in neutrino propagation in matter. 

Recently, this problem of searching for higher dimensional operators without 
charged lepton constraints has come to conclusion; It has been proved that the above 
two possibilities are unique in dimension six and eight operators, respectively, if one 
wants to avoid the charged lepton constraints at the tree level 

Now, the question is: Is avoiding the tree level constraint sufficient to be free from 
the charged lepton constraints? It is known that the answer is N0^22l_ Namely, the 



dressing by W and Z bosons can produce four charged lepton processes which leads 
to highly restrictive bound in some channels, in particular on \el'^^\. See (new 
version) for details of the type of bound, and for a summary of the other constraints 
on NSI. 

It turned out, however, that it was NOT the end of the story. Biggio, Blennow 
and Fernandez-Martinez have recently pointed out that the bounds on NSI have 
to be relaxed to a large extent, a factor of ~ 10^. Because of the antisymmetric 
nature of the dimension six operator ([2]), the contributions of diagrams with different 
flavor indices tend to cancel and add up to zero in the limit of neglecting the lepton 
massesE Turning on lepton masses leaves the contribution of the order of {rrii/Mw)'^ ■ 
Notice that it is the unique dimension six operator which is free from the tree-level 
four charged lepton counter part so that we have to live with it. Another significant 
feature is that the bound on goes away, or in other word, it must vanish by 
construction of the operator ([2]). 

I would like to note here that an important feature is hidden behind my shallow 
description of their results; SU(2) gauge invariance. That is, imposing the gauge 
invariance is essential to obtain gauge invariant results of lo gari thmically divergent 
terms of the one-loop diagrams. Go to the original reference for more complete 
understanding, in particular on the meaning of quadratically divergent terms. 



4. NSI in Neutrino Propagation in Matter 



In the rest of my presentation I discuss NSI effects in neutrino propagation in mat- 
ter. It should be remarked, however, that NSI effects are present also in production 
and detection processes of neutrinos, so that my discussion is obviously incomplete. 

To summarize its effects on neutrino propagation it is customary to introduce the 
e parameters, which are defined as Ea/s = Hf,p '^^{S^ where rif is the number density 
of the fermion species / in matter. Notice that only the vector combination of the 
NSI can be probed when we discuss neutrino propagation in matter. Approximately, 



the relation Ea/B — J2p + 3 + 3 e^^j holds because of a factor of ~3 larger 
number of u and d quarks than electrons in iso-singlet matter. Notice that with the 
dimension six operator ([2]) part of the first term, the ones from and is absent. 
I, however, choose to proceed with a generic framework. 

Using the e parameters the neutrino evolution equation which governs the neutrino 
propagation in matter is given as 
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Here, I restrict ray discussion into the dimension six operators. Bu t, it is fair to note that this 
discussion is much more relevant for the dimension eight operators. See 1241 _ 



where U is the MNS matrix and a = 2\/2G pneE ^ where E is the neutrino 
energy and Ue denotes the electron number density along the neutrino trajectory in 
the earth. Amf^ = mf — with neutrino mass (i = 1 — 3). Notice that the phase 

of e parameters may provide new source of CP violation 1^^. Another important point 
is that complexity of the system in^^ would lead to confusion in determination of 
the mixing and the NSI parameters 1^^^'. 



5. Perturbation Theory of Neutrino Oscillation with NSI 

Obviously, I am a newcomer to the field of NSI. When I started to work on this 
topics I tried to understand the features of neutrino oscillations with NSI. Alas, I 
found that not so many things are known. The questions I would like to know the 
answer were: 

• From the experience in neutrino oscillation with standard interaction (SI) I 
would expect that the appearance channel z/g — > z/^ (or, Vt) has great 
sensitivities to tiny effects of NSI. Then, the natural question is; Which NSI 
elements of Sap in (EI) give the dominant contribution to P{ve ^^i)? Or, more 
concretely, how large is the contribution of e.g., Eefj., £er, and to P(z^e 

• What about the disappearance channels though they may be less attractive? 
Namely, what is the size of contributions of Eee in P{i'e —>■ J^e)? What is the 
relative importance of Ee/j,, Eer, ^^n, and E^r in -P(t'^ t'lJ? 

I was amazed by the fact that apparently nobody knew the answers to these 
questions!^ I believe that the questions are not only due to academic interests. It is 
because I think that treating the full system ([6]) is really necessary. Though people 
(including myself) do make approximations of ignoring some elements keeping only 
a few of them, but they do so without good reasons. It is even more so now because 
most of the stringent bounds on Eajj based on lepton processes, the model-independent 
ones, are gone. Only when we recognize the correct theory at high energy scale we 
can be sure that the approximation he/she is making is the correct one. 



5.1. t Perturbation theory 

To answer these questions and to have a global bird-eye view of neutrino oscillation 
with NSI we have formulated a perturbative framework^l^. Unfortunately, there is no 
unique framework because we still do not know the value of ^13, though the bound on 

It might be a too strong statement, given the fact that so many people are working in this field. 
Any comments are welcome. In fact, it appears that answer to these questions were known at least 
partly by Jacobo Lopez-Pavon in UAM, Madrid, though the result was unpublished. 



it exists"^. The only parameter which we know to be small is the ratio — 0.03. 
Therefore, we take an ansatz 



e=- — 2 Si3~£q/3~10 {a,(3 = e,n,T) (7) 



to formulate our perturbation theory, which we called the "e perturbation theory" 
in [13. In doing so I assume . is of order unity, anticipating very long-baseline 



neutrino experiments such as neutrino factory or the beta beam We do 
not take (^-^ — S23) as an expansion parameter because a rather large range is cur- 
rently allowed and the situation will not be changed even with the next generation 
experiments . 

The e perturbation theory based on ([7]) is a natural generalization of the frame- 
work taken by Cervera et al. ^ for systems with SI into the one with NSI. Note 
that the Cervera et al. formula (which we call the SI second-order formula) is the 
most widely used perturbative formula to discuss various aspects of neutrino oscilla- 
tions. We derive the NSI second-order formula which generalizes the SI second-order 
formula to the case with NSI to have an overview of the neutrino oscillation phenom- 
ena in systems with NSI. A different but related perturbative approac h to neutrino 
oscillation with NSI has been studied in references 1 32|33 | 34 | 35 | 17 | 36|37 | 

In passing I have a few remarks on ^13. It is a big question whether 9^-^ falls into 
the range which can be explored by the next generation accelerator '^^ ' and the 
reactor 1^ experiments. In fact, I argue sometimes rather strongly that 61^ must be 
large, for example in^^jf] The belief is one of the motivations for my works which 
proposed reactor measurement of ^13 and superbeam measurement of lepton CP 
violation!^. Nevertheless, I am a pessimist here with the ansatz ([7]). Well, the 
reason why I take the ansatz of small ^13 is that it is the only natural perturbative 
framework of neutrino oscillation. For instance, the appearance oscillation probability 
P{ve — i^^i) consists only of order terms. If I take a different ansatz S13 ~ ^/e = 



(which roughly correspond to the Chooz limit 123)^ the terms in P(z/e u^) 

3 

do not scale uniformly and we would have to keep terms of order ea to include effects 
of CP violation. It necessitates to keep order 5^3 terms. 



5.2. NSI second-order formula; ly^-related sector 

How can one go to the NSI second-order formula from the SI second-order formula? 
Though the task might look formidable, it is in fact trivial in i/g-related channels! 
What is necessary is to make replacements in the atmospheric and the solar variables 

° The basic reasoning for my belief is simple:'^^' The MNS matrix is the product of the two 
unitary matrices which diagonaUze the neutrino and the charged lepton mass matrices. The two 
angles in the MNS matrix are known to be large. Then, why should the third one extremely small? 



in the SI second-order formula and that's it: 



Sl3 Si3 + is23eef^ + Caster 6*', 

a a 

C12S12 C12S12 \- C23£eiJ. — S23£eT- (8) 

a a 

It is very easy to understand why the particular combinations of e parameters come 
in to the atmospheric and the solar variables, respectively. It is well known that 
in doing perturbation theory the convenient basis is the tilde basis H = U2sHU23- 
The combinations of the e parameters are the ones that appear in the NSI part of 
He3 and He2, analogues of "sin6'i3" and "sin6'i2". (See equation (15) in^^.) 



The resultant NSI second-order formula in channel reads 




X . / Am|i \ . /aL\ ( .AmlL 

P[^e ^fi) = 4 C23 ( C12S12— ^ h C23£eM " S23£er ) Sm ( — 1 CXp I -I 

+ S23 [Si^e — h S23£ef, + C23£^eT j (^ ^^2 _^ J 

1 hope the readers are convinced of my claim that the formula is surprisingly simple 
in its form. P(z/e — ^ Ut) can be obtained by doing the transformation C23 —>■ —S23 
and S23 —>■ C23 in P(z/e u^), but undoing any transformation in the generalized 
atmospheric and the solar variables defined in ([8]). 

A notable feature in ([9]) is that only the elements Eefi and Eer appears in the NSI 
second order probability formula. Because of the decoupling of the other e's it is in 
principle possible to determine and Eer together with the SI parameters 6*13 and 6, 
6 real parameters including phases. If one carries out this task by rate only analysis 
we need measurement of the oscillation probabilities in the following three channels 
i^e i^^i, ^ i^T, and z/^ — > Uf. and their antineutrino counterpart. 



5.3. NSI second- order formula; — v^. sector 

In i'fj_ — Vr sector the situation is different. The NSI dependent piece in the 
oscillation probabilities P{y^ — >■ i^r), P{y^l J^fi), and P(^'r — >■ i^t) is universal. 
See for explicit expressions. Because of this feature one cannot determine all the 
relevant NSI elements E^r and e^^ — Err (only the difference can be measured), 3 
unknowns, by the rate only analysis. 

The reasons for such curious feature of the universal NSI dependent term is simple 
to understand. By unitarity it follows that 

P(Ur^Ur)+PiUr^U^) = l-P^u^^U,). (10) 



We note that Pi^i'^ Ve) and P(z^r ^e) do not contain e^^, Err, and to 
second order in e. Then, it follows from the first equation in ffTOl) that ^ 
i^T',£fj,ii,£fiT,£TT) = —P{t^fj, i^ti'T^tiiiy^fiTT^rr)- Notlclng that the terms related to 
e's in the — Ur sector are T-invariantij the relations P(z/t- — >• lyr] ^t^-fj-^ ^^J■r^ ^tt) = 
— P(f^ Pr] £^lfJ., £^iT, £tt) must also hold. Therefore, the Eap (a, P = fi,T) dependent 
term in the three channels are all equal up to the over-all sign. The necessity of 
spectrum analysis is obvious to determine all the NSI and the SI parameters. 



5.4- Summary table 

To answer to the questions raised above I present below the summary table. 
One of the features in Table. [1] which requires comment is that in the last column. 
In standard neutrino oscillation only with SI the matter effect comes in into the 
oscillation probability only at the second order in e, the property dubbed "matter 
hesitation" in^I^. It is the reason why it is so difficult to detect the matter effects in 
many accelerator neutrino experiments including NOz/A ^SHl ^pj^g matter hesitation 
is a highly nontrivial property because we treat the coefficient a (I mean, a/Arn^i) 
as of order unity. For example, there exists first order a dependent term in the S 
matrix, but it does not survive in P(z/e ^e) because it enters as a phase factor. 
Notice, however, that its validity relies on the particular framework of perturbation 
theory. For a (simple!) proof of this property see^^. 



Table 1: Presented are the order in e (~ 10"^) at which each type of Sap {oi,/3 = e,/i,T) and 
a dependence (a is Wolfenstein's matter effect coefficient) starts to come in into the expression of 
the oscillation probability in e perturbation theory. The last column is for the a dependence in the 
standard oscillation without NSI. The order of e indicated in parentheses implies the one for the 
maximal 623 in which cancellation takes place in the leading order. 



Channel 


See £e/i ^er ^/it ^/J-fJ- ^tt 


a dep. (NSI) a dep.(SI) 


a = e, ii,T 


g3 g2 g2 g3 g3 g3 






e3 e2 ^2 ^1 gl^g2) ^1(^2) 


el £2 



One of the implication of the matter hesitation is that See comes into the oscillation 
probability at order in all channels, as indicated in Table. [U It is because See is 
nothing but a small shift of the matter effect coefficient a. Because of this property it 
is very difficult to measure e^e in long-baseline experiments. It should be remarked, 
however, that assuming that the other NSI elements are vanishi ngly s mall it can be 
measured in a great precision of a few % at a neutrino factory But, I must 



We emphasize that this feature itself is highly nontrivial, and can be realized only by an explicit 
computation. 



note that it is only true under the assumption that the earth matter density along 
the neutrino trajectory is accurately known. 

Another notable feature in Table. [1] is that there exists first order term of NSI 
element e^,- (and e^^ — Srr if 6*23 7^ f) in the z/^ — Vr sector. Clearly, they are due to 
direct transition caused by these NSI elements. In fact, rather high sensitivities for 
determining and e^^ — e^r observed in atmospheric ^Sl a^ci future accelerator 
neutrino analyses are benefited by this feature. 



5.5. SI— NSI confusion 

The structure of the oscillation probability ([9]) in which e^a enters into the ex- 
pression only through the generalized solar and the atmospheric variables ([H]) is one 
of the most significant features of the oscillation probability with NSI. It also implies 
that, when the parameter determination is attempted, there exists severe confusion 
between the SI and the NSI parameters. Hence, our result may be regarded as an 
analytic proof of the general "confusion theorem" in e perturbation theory. 

The uncovered structure may be helpful to formulate a strategy of resolving the 
confusion, because ([8]) clearly dictates which SI parameters will be confused by which 
NSI variables by which way. Notice that our confusion theorem is quite different in 
nature from the one proved which is confused with the NSI elements in 

production and detection processes. 



5.6. Parameter degeneracy in neutrino oscillation with NSI 

It is now well understood that phenomenon of parameter degeneracy, existence of 
the multiple solutions, occurs in neutrino oscillation measurement of SI parameters 
|51 52 53) Because of the large number of unknown (i.e., to be determined) param- 
eters (2 standard and 8 NSI ones) the parameter degeneracy in the full system is a 
formidable problem to work out, even under the approximation of ignoring NSI effects 
in production and detection of neutrinos. 

Yet, it was possible to recognize a completely new type of degeneracy . Let us 
denote the generalized atmospheric and the solar variables in ([8]) in section 5.2 as G-t 
and Se"*'', respectively. Then, if there is a solution \Q^±\ and jS*^^-*!, then the second 
solution lel^^l = and = ^I±\Q^l\ exists. (See for definitions 

of X± etc.) It can be called the "atmospheric— solar variables exchange" degeneracy, 
which arises because of large number of unknown parameters in the solar and the 
atmospheric variables. Of course, it does not survive when NSI is turned off because 
there is no solar degrees of freedom (as to be determined parameters) as can be seen 
in dH]). 

What is the right way in this difficult problem of degeneracy in systems with 
NSI? As a first step, we have worked out the problem in a region where the matter 



effect can be tr eated as a perturbation. For early references of matter perturbation 
tfieory, see It is known that analysis of t he parameter degeneracy becomes 

particularly transparent in this setting l 52 | 56 | 57 [5Hli 

Our analysis is most transparent in the "discrete" degeneracy, the sign-Am|^ and 
the octant ones. Let me describe first the sign-Amg^ degeneracy. One can show that 
the NSI dependent terms in the oscillation probability P(z/e — ^ i^a) {(^ = fj,, r) to first 
order in a is invariant under the transformation 

6 71 — 6, 

4>ea 27r-0ea, (H) 

while keeping 613 and l^eal fixed. It nicely complements the discussion in^^ and it 
indicates that there exists a new (approximate) solution with differing sign of Am'^i. 
It is worth to note that the invariance is true only if the CP phase of NSI element is 
involved in the transformation ( ITTi) . 

Similarly, one can show that there is another invariance under the transformation 
(assuming 6^23 7^ f) 

C23 ~^ S23, 
S23 ~^ C23, 

It means that the ^23 octant degeneracy prevails in the presence of NSI, and actually in 
an extended form which involves NSI parameter e^^ — Err- Since this NSI parameter 
decouples from P^i'^ I'e) to second-order in e, the presence of the 6*23 octant 
degeneracy remains intact when the NSI is included though values of the degenerate 
solutions themselves are affected by the presence of Eea- Thus, we have shown that 
the parameter degeneracy survives the presence of NSI provided that NSI elements 
are "actively involved" in the degenerate solutions. 

The other salient feature of degeneracy in the mass perturbative regime is the 



iiroi 

property called the "decoupling between degeneracies" . We have revisited this 
issue in the systems with and without NSI. The conclusion obtained in is that 
the decoupling between the sign-Amg^ and the 6*23 octant degeneracies holds with 
and without NSI. On the other hand, decoupling between them and the intrinsic 
degeneracy does not holds with and without NSI, partly correcting the conclusion in 



The only exception is the special setting at the oscillation maximum, or more 



precisely, the shrunk ellipse limit 1^^, for which the decoupling holds. 



5. 7. Parameter degeneracy; An example 




Figure 1: Allowed regions in the 0e^t — \sefj\ plane (left column), (p^^ — 5 plane (middle column) and 
sin^ 2^13 — 5 plane (right column) corresponding to 2 and 3 cr CL obtained for the system with e^,^ . 
Panels in the upper 2 rows (3rd row) correspond to the case where only a single detector at 3000 
km (7000 km) is taken into account, whereas the ones in the 4th row correspond to the case where 
results from the two detectors are combined. The input parameters are taken as: sin^ 2^13 — 10~^, 
5 = 0, Ice^l = 5 X 10"'^ and 0e/j = 57r/4 (indicated by the green asterisk), and the mass hierarchy 
is normal. For the case where only the single detector at 3000 km is taken into account, allowed 
regions exist not only for the normal mass hierarchy regime but also for the inverted one, as shown 
in the panels in the second row. 



An example of the parameter degeneracy is presented in Fig. [TJ This is one of the 
examples found in doing the work but is not presented in the reference. It clearly 
demonstrates the existence if th e Am^ i-sign flipped and the intrinsic degeneracies, 
the natural extension of the one 1 ^^ ^ ^^^ to the systems with NSI. As pointed out in 
section 5.6 the phase (pe^i of the NSI element is indeed heavily involved though 
the relation fill I) which is valid in the mass perturbative regime does not quite hold. 
In this example (as well as in the one which is presented in^^SJ)^ f^j- detector 
measurement successfully lift the Amg^^-sign flipped degeneracy, but not and the 
intrinsic one. 

How robust is the degeneracy in system with NSI? It is a difficult question to an- 
swer in general. But, it is worth to remark that sometimes the degeneracy is extremely 
hard to solve because the energy spectra corresponding the degenerate solutions are 
so similarll This fact is demonstrated in Fig. [2] in the e^t system with the particular 



8 It is often the case that in systems only with SI that the intrinsic degeneracy can be "easily' 




Figure 2: Energy spectra of the oscillation probability for a system with non-zero NSI element Ser 
corresponding to the same Am§]^-sign intrinsic degeneracy (left panel) and the flipped Am§]^-sign 
degeneracy (right panel). 



values of the SI and NSI parameters. See Fig. 13 in^lSi for the corresponding figures 
for the system with exactly the same feature. 

6. Discriminating CP violation due to SI and NSI phases 

I emphasize that one of the most important features of the system with NSI is 
the coexistence of two kind of CP violation , the one due to 5 in the MNS matrix 

the lepton ic vers ion of the celebrated Kobayashi-Maskawa (KM) phase in the 
CKM matrix 1 ^"^ ! ^^^ for quarks, and the other which come from the phases of NSI 
elements. Knowing the nature of CP violation seen in any kind of experiments is 
of decisive importance because of many reasons, in particular for possible connection 
to leptogenesis scenario currently the most promising one for baryon number 
generation in the universe. 



6.1. Two -detector setting 

I have discu ssed in the last year in Venice the possibihty of resolving the ^13— NSI 
confusion I^SI^ by the near (3000 km)— far (7000 km) two-detector setting in neutrino 
factory based on our work \^ . I use the same setting to examine the question of 
whether the two-phase confusion can be resolved For a related work on the same 
subject see . The similar question of distinguishing two kind of CP violation in 

resolved by the spectrum analysis. See for example, in the case of T2K or T2KK settings l 57 | 58 |i^ 



the context of "unitarity violation" approach has also been investigated inESEU 
Setting of the second detector at around the magic baseline ~ 7000 km was moti- 
vated by high sensitivity to the matter effect . It is in concordant with the similar 
two detector setting in a neutrino factory as a degeneracy solver 1 51 1 6716 8} _ rj.^^_ 
detector setting has been proposed in neutrino experiments in a variety of contexts^ 
The basic idea in the present context is to seek complementary role played by the far 
detector. 



6.2. Use of the bi-probability plot 

Though I told you that derivation of the NSI second-order formula of the os- 
cillation probability is simple, it means neither that the dynamics of the system is 
simple, nor it is easy to understand. What makes the system so complicated is the 
very existence of two CP violating phases, 6 and the phase 0ea of the NSI element 
Eea = keole*'^'"" (« = /U, t). In my talk, therefore, I report the work done with only a 
single type of NSI, either first step of understanding the features of 

neutrino oscillation with NSI. 

How complicated is the system with NSI? Seeing is believing. Presented in Fig. [3] 
is the bi-probability plot in P(z/e z/^J — P(z/e ^^i) space - ■ but by varying the 
two phases, 5 and 0e^i- As you see the ellipses move around in the plane such that 
the whole triangular region is (almost) swept over. So "anything can happen" with 
the two phases. We have characterized the behavior of ellipses as rotating ellipses 
in ; An ellipse drawn by varying the phase A rotates when the other phase B is 
varied. Because of the behavior of the probabilities rich phenomena such as confusion 
between SI and NSI parameters and the parameter degeneracy are expected. 

In fact, various viewpoints have to be involved to really understand features of 
neutrino oscillations with NSI and the sensitivities to the NSI elements, \eea\ and (pea 
(a = /i, t), and the SI parameters, 5 and sin^ 2^13, to be achieved by the detectors at 
L = 3000 km and L = 7000 km separately and in combination. They include: 

• How prominent is the synergy between the near and the far detectors for deter- 
mination of SI and NSI parameters? How it differs between the systems with 

^ Here is a brief history of two-detector settings in contemporary neutrino experiments: It was 
proposed as an appropriate setting for measuring CP violation in the context o f lo w energy 
superbeam experiment ' . In a quite different context of reactor measurement of 613 the two- 
detector setting is the standard one to guarantee the near-far cancellation of systematic errors. It 
has triggered interests in the world wide scale , and led to the several international collaboration 
experiments SU . The idea has also been applied to the Kamioka-Korea two detector complex with 
an upgra ded ne utrino beam from J-PARC to determine the mass hierarchy as well as discovering CP 
violation I 57 | 58 [ _ "unifies" the two aspects of near-far cancellation and syn ergy between the two 
detectors, and can ser ve for a possible upgrade option of the T2K exp eriment . For an overview 
of T2KK, see e.g.,1^, and for a review of the two-detector setup ^^^F, 




Figure 3: Bi-probability plots in P{i'e v^) — P{ve v^) space dX L — 3000 km, for E ^ 20 GeV, 
with sin'^ 26*13 ~ lO^'^ and Eefj. = 5 x 10"'^ computed numerically using the constant matter density 
p ~ 3.6 g/cm"^ assuming the electron number density per nucleon of 0.5. The both axes is labeled 
in units of 10^*. The values of the parameters taken are £e/i sin^ 2^13 



EepL and ^er? What about dependence on values of the parameters, in particular 
on the size of NSI elements? 

• How can the two-phase confusion be resolved? Does the answer depend on 
which NSI element is turned on? 

• What is the nature of the parameter degeneracy in system with NSI, and 
whether it can be resolved by the two-detector setting? If so how can it be 
realized? 

These points are fully discussed in our paper ^^3. 

Here, we make comments only on a puzzle. The difference in sensitivity to NSI 
has non-trivial features. At relatively large values of NSI, e^^ = 10~^ and Ser = 10~^, 
the size of the ellipses are similar in size. But, the sensitivity is in fact very different 
between the two systems; The one in the Eefi system is much higher than that in the 

system. On the other hand, the parameter degeneracy is much severer at the near 
detector in systems with ^e/x compared to the ones with shown in the tables. 

A unified understanding of the puzzling features becomes possible once one draws 
the bi-probability plots by varying S for several different values of and ^13, as done 
in Fig. m The degeneracy is severer in the ^e^i system because of the more dynamic 
behavior of the bi-probability ellipses as shown in the left panel of Fig. HI Because 
the ellipses can locate themselves essentially everywhere in the bi-probability space 
there are chances that fake solutions can be produced at points far apart from the 
true solution. What about the difference in the sensitivities? We observe in the right 



CP elipses for fixed if) and varying 5 



CP elipses for fixed ^ and varying 5 




P(v^^v^ ) [xlO"'] P(v^^v^ ) [xlO"'] 



Figure 4: Bi-probability plots drawn by continuously varying S for four different values of (f). 

panel of Fig. Hlthat the ellipses in the Eer system remain in much more compact region 
when 6 and 613 are varied. Because of the finite resolution of the experimental data 
it appears that the dense concentration of the ellipses with different parameters leads 
to merging of many degenerate solutions. It probably explains lack of the sensitivities 
and at the same time much less frequent degenerate solutions in the Ser system. 

7. Discovery Potentials 

This is the appropriate point to discuss the discovery potential of various quanti- 
ties, l^e/^l, kerl, 4>e^i, 0er, as wcU as the Standard parameters S (and 6*13 in principle) 
and the neutrino mass hierarchy. In this report we focus on CP violation caused by 
the lepton KM phase 6, and the pha se of NSI. The discovery potential for the rest 
of the quantities are discussed in^^^. 

By the way, I remind you that all the figures presented in this manuscript are 
new, i.e., no single figure which is identical to the one in^^^. To keep this tradition I 
will always give in this manuscript the sensitivity regions calculated with the inverted 
mass hierarchy as input. Notice that all the sensitivity plots given in^^ are calculated 
by taking the normal hierarchy as input. Great thanks to Hiroshi Nunokawa for his 
efforts to prepare them for this manuscript. 

While I do not give any details of the quantitative analysis in this manuscript 

(for which see ) , it should be remarked that all the systematic errors as well as 

backgrounds are ignored in our analysis. Here, I explain the reasons for this choice in 

my own language. The dimension six operator that ca n g ive neutrino's NSI without 
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producing unwanted four charged lepton NSI is unique , the anti-symmetrized one 
given in Eq. ([2]) in section 3. Then, unless someone is able to show that only the 



dimension six operator naturally arises in a certain class of models of new physics at 
TeV scalefl we must prepare for search for the dimension eight (or higher) operators, 
with the size |£:| ~ 10~^ assuming no extra suppression^] If it turned out to be the case 
one must think of the experimental technology which can accommodate this request. 
Since its realization is not known, we invented a model experiment using the neutrino 
factory setting with an ideal detector for which the systematic errors and background 
are ignored. 



7.1. CP violation due to NSI 



Let us start with the sensitivity to CP violation caused by the phase of NSI. In 
Fig. [5] and Fig. [6] the regions sensitive to non-standard CP violation due to NSI are 
presented in — |£:| space. I these regions one can detect non-standard CP violation 
by NSI (0 7^ and 7^ vr) at 2(T (red thin lines) and 3cr (blue thick lines) CL. Fig. [5] 
and Fig. [6] are for the and the e^r systems, respectively. 




Figure 5: Regions where the non-standard CP violation caused by (j)e^ ^ or (t)e^ ^ tt can be 
estabUshed for the case sin^ 26*13 = 10~^, (5 = tt (left panel) and 5 = 37r/2 (right panel). The 
inverted mass hierarchy is assumed as the input. 

By comparing the Fig. [5] and Fig. [6] to Figs. 18 and 20 in^^^, respectively, one 
notices several notable differences between the normal and the inverted mass hier- 
archies. In the Eefj^ system the sensitivities to non-standard CP violation at 5 = vr 
are significantly worse both at the near detector (3000 km) and the near-far (7000 
km) combined in comparison to those obtained with the input normal hierarchy. The 
results at 5 = 37r/2, however, are very similar to the case of normal mass hierarchy. 

' From the reasoning below, I think it important to pursue this possibility. 

j I have assumed throughout this report that discussions on NSI in the lepton sector applies to 
the operators which involve neutrinos and quarks. Though I think it reasonable it can be subject to 
criticism. 
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Figure 6: The same as in Fig.[5]but for the non-standard CP violation caused by 4>eT or (per tt. 

A somewhat curious behavior seen in the upper-right panel of Fig. [5], no sensitiv- 
ity to non-standard CP violation at the maximal CP violating inputs, = 7r/2 and 
= 37r/2 is explained as a consequence of the parameter degeneracy, the one called 
the (p degeneracy in^I^. 

In the Eer system the sensitivities to non-standard CP violation are similar to the 
normal hierarchy case. The most notable difference is in the 6 = it case; At the near 
detector the sensitivity to non-standard CP violation is a bit worse than that of the 
normal hierarchy case, but curiously enough it is a little better when the far detector 
is combined. It can well be the case because the features of synergy between the two 
detectors are highly nontrivial ^^^^^ . 



7.2. Standard CP violation 

Let us go back to the sensitivity to CP violation due to the KM phase 6. In Fig. [7] 
the regions sensitive to the standard CP violation due to 6 are presented in the system 
without NSI. They are significantly worse than the normal hierarchy case given in 
Fig. 22 in . Most probably, it is due to relatively smaller number of events in the 
antineutrino channel. Notably a peninsula like region without sensitivity develops 
from S ~ O.Stt to (5 ~ O.Svr. 

The gross features of the sensitivity regions remain unchanged even when the NSI 
degrees of freedom is turned on, as can be seen in Fig. [HI The sensitivities to the 
standard CP violation are slightly worse compared to the normal hierarchy case given 
in Fig. 23 intlSl. 

I must warn the readers that the sensitivity contours are unstable to inclusion of 
the systematic errors and backgrounds. Yet, I suspect that these sensitivity regions 
are similar to the ones obtained with the apparatus which can explore the NSI in the 



whole region down to |£:| 
the future NSI search. 



10 ^, the sensitivity which I argued to be necessary in 
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Figure 7: Sensitivity to discovery of standard CP violation. Here no effect of NSI is assumed in 
the input data or in the fit. The upper panel shows the case where only the detector at 3000 km 
is considered, whereas the lower panel is the case corresponding to the combination of detectors at 
two different baselines. The inverted mass hierarchy is assumed as input. 
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Figure 8: The similar plots as in Fig. [7] but with non-zero NSI allowed in the fit; The input data 
was generated assuming the inverted mass hierarchy without NSI but non-zero values of ee/i (left 
panel) and Ser (right panel) were allowed in the fit. 



8. Conclusion 

After reviewing the theoretical progresses on NSI recently made in section 3, I 
tried to explain the works done by our group on theoretical and phenomenological 
aspects of hunting the NSI in sections 5-7. I guess the former contributed to illuminate 
the global features of neutrino oscillation with NSI including method for parameter 
determination and recognition of the parameter degeneracy with NSI. While in the 



latter we have investigated the problem of discriminating the two kind of CP violation, 
one due to the standard KM phase and the other by phases of NSI elements. 

It appears that the near (3000 km) — far (7000 km) two detector setting in 
neutrino factory does have a rather high sensitivity to explore \Se^\ to ~ 10~^ in a 
lucky region of and to ~ several xlO""^ in most of the region of 0. The sensitivity 
to leerl is lower but still it can be explored to ~ 10^^. See Figs. 14-17 in^I^. The 
sensitivity to CP violation is also very good; The one due to NSI phase can be probed 
to l^e^l to ~ several xlO~^ in most of the region of 5 and (/)e^. They are close but 
not quite the lower end of the required full region for exploration of NSI due to 
TeV scale new physics. Moreover, our estimation ignores the systematic errors and 
backgrounds, and hence is overly optimistic one. 

How can this situation be overcome? Honestly, I don't know the answer. However, 
a few comments may be made: 

• Effects of NSI in production and detection of neutrinos, which are completely 
ignored in our works, can be of great help. 

• One can formulate the good enough reasoning to convince people that search 
for neutrino's NSI to \e\ ~ 10~^ is sufficiently informative to signal new physics 
at TeV scale. 

However, I guess the former possibility is not easy to be realized. Nonetheless, I 
would like to recall that no discovery done in the past was an easy one. Also, upon 
identification or grasp of the new physics we should obtain the clearer view. A second 
eel can be a big one! 

9. Bibliographical Note 

Given the large number of references devoted to the subject of NSI it is not easy 
to find the appropriate one, in particular for a newcomer as I was sometime ago. 
Therefore, I tried some efforts to collect them here with classification by subjects. 
Yet, it is extremely difficult to find all of them, and therefore, the list should be 
considered as an incomplete one. I would like to apologize to those who will find 
their references missed. The categories I use are: 

• Accelerator neutrinos; neutrino factory l ^^ l ^^ ^ '^^ l ^'^ l ^^ l '^'^K 

• Accelerator neutrinos; excluding neutrino factory I B5 | 50 | 74 | 

• Atmospheric neutrinos ^I^^^^^^. 

• Reactor or spallation source neutrinos or low-energy scattering . 

• Solar neutrinos . 



Astrophysical neutrinos 
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